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Motivated by the studies of the superconducting pairing states in the iron-based superconductors, we analyze 
the effects of Brillouin zone folding procedure from a space group symmetry perspective for a general class of 
materials with the PA/nmm space group. The Brillouin zone folding amounts to working with an effective 
one-Fe unit cell, instead of the crystallographic two-Fe unit cell. We show that the folding procedure can be 
justified by the validity of a glide reflection symmetry throughout the crystallographic Brillouin zone and by the 
existence of a minimal double degeneracy along the edges of the latter. We also demonstrate how the folding 
procedure fails when a local spin-orbit coupling is included although the latter does not break any of the space 
group symmetries of the bare Hamiltonian. In light of these general symmetry considerations, we further discuss 
the implications of the glide reflection symmetry for the superconducting pairing in an effective multi-orbital 
t — Ji — J 2 model. We find that the PA/nmm space group symmetry allows only pairing states with even 
parity under the glide reflection and zero total momentum. 


I. INTRODUCTION 

The iron-based superconductors form a large family of ma¬ 
terials which exhibits considerable diversity in their lattice 
structures. Examples include the 1111, 111, and 122 iron 
pnictides, as well as the 11 iron chalcogenides. The struc¬ 
tural unit cell of most of these superconductors consists of 
two Fe and two pnictogen/chalcogen atoms and it is typically 
labeled as a 2-Fe unit cell. In momentum space, the corre¬ 
sponding physical Brillouin zone (BZ) is usually referred to 
as the folded BZ (FBZ). However, many theoretical studies 
have been based on tight-binding models defined on an effec¬ 
tive ’’unfolded BZ” (UBZ) of 1-Fe unit cell with an implicit 
equivalence between the former and the FBZ. The mapping 
between the effective 1-Fe UBZ and the physical 2-Fe FBZ is 
via a BZ folding procedure. In view of the widespread use of 
this mapping we believe it is important to better understand 
the necessary conditions for its employment. 

An early discussion of such a BZ folding procedure for 
1111 iron-pnictide superconductors was provided by Fee and 
Wen [1], They noticed that a single Fe-As plane contains a 
glide reflection symmetry which consists of a fractional unit 
cell translation followed by a reflection about the Fe-plane 
(see below). This symmetry can be used to define a pseudo¬ 
crystal momentum which can label the single-particle wave- 
functions with different parities under the operation. The 
last step in turn allows the definition on an unfolding pro¬ 
cedure from the 2-Fe FBZ of the physical momentum space 
to the 1-Fe UBZ of the pseudo-crystal momentum space. 
More recently, a thorough group theoretical description of the 
electronic structure in the iron-based superconductors with a 
PA/nmm space group symmetry has been given in Ref. 2. 
There, it was found that the glide-reflection symmetry clas¬ 
sifies the Bloch states near the Fermi level and puts strong 
constraints on the low-energy effective model of the system. 
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The immediate motivation for the current work has come 
from the strong-coupling approach to superconductivity in the 
multi-orbital models for the iron-based materials. Supercon¬ 
ducting paring in this approach has been studied by using 
t, — Ji — J -2 models with three or more 3 d orbitals, involv¬ 
ing at least the 3d xz , 3 d yz and 3 d xy set, in an effective 1-Fe 
unit cell [3-5]; such studies have been motivated by both the 
multi-orbital nature of the electronic structure and the bad- 
metal behavior of the parent compounds [4], As we will dis¬ 
cuss, the Bloch states formed directly from 3d xz and 3 d yz 
Wannier states are odd under the reflection about the Fe plane, 
whereas those associated with the 3 d xy orbital are even under 
this symmetry operation. In the former case, the unfolding 
procedure becomes trivial as both types of Bloch states map 
onto the same quasi-crystal momentum in an unfolded BZ. 
By contrast, states derived directly from the 3d xy orbital are 
mapped onto a quasi-crystal momentum which is shifted w.r.t. 
the other two orbital states. Given this, the validity of calcu¬ 
lations done directly in a 1-Fe unit cell comes into question 
and a careful examination is required to establish whether a 
more involved 2-Fe unit cell basis should be used instead. It 
appears that this potential issue is not restricted to the particu¬ 
lar t — J\ — J2 case. Indeed, within a more general but related 
context, different ways of taking into account the two inequiv¬ 
alent pnictogen/chalcogen atoms have recently been proposed 
[6, 7]. A consensus does not appear to be reached since even 
the lattice symmetries considered in these studies are quite dif¬ 
ferent: Ref. 6 considered a C,\ v point-group symmetry about 
each As atom, while in Ref. 7 a local £>4 symmetry about each 
Fe atom was taken into account. The proposed effective tight- 
binding models and the superconducting pairing symmetries 
in these two works are also quite distinct. Indeed, while Ref. 6 
proposed a spin singlet, orbital triplet d + id A \ g pairing (de¬ 
noted as “TAO pairing”) as a consequence of the two inequiva¬ 
lent As atoms. Refs. 7-9 discussed an s-wave odd-parity pair¬ 
ing with nonzero total momentum, labeled 77-pairing. These 
proposals not only question the validity of using the an effec¬ 
tive 1-Fe UBZ, they also point to the need for a clear link be¬ 
tween the superconducting pairing and lattice (space group) 
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symmetries. This is particularly the case for the 77 -pairing, 
where it is still uncertain whether such pairing is allowed by 
the space group symmetry [ 8 - 11 ]. 

We believe that part of the confusion regarding the pair¬ 
ing symmetry is related to the issue we mentioned in the be¬ 
ginning of our paper: To what extent can a model for iron- 
based superconductors defined on the 1-Fe UBZ reproduce 
results consistent with the one defined on a 2-Fe FBZ? Here 
we study this problem for systems with a P4/nmm space 
group (which include the 111 and 1111 iron pnictides and the 
11 iron chalcogenides). The stringent conditions imposed by 
the space group constitute our starting points in the analysis of 
the validity of the 1-Fe unit cell formulation. To our knowl¬ 
edge there have been only a few attempts [1, 2, 12, 13] at 
placing this procedure on a firmer footing. More specifically, 
we wish to give a better description of the notion of the glide 
reflection symmetry within a more formal group-theoretical 
context and provide some justification for it’s implicit use as a 
general symmetry which can be used to label Bloch states of 
arbitrary momenta. 

The glide reflection operation is a part of the P4/nmm 
space group, and is characteristic of the non-symmorphic na¬ 
ture of this group. As we will see, it plays an important role 
in establishing, in the absence of any spin-orbit coupling, the 
validity of using the 1-Fe UBZ for both the single-particle dis¬ 
persion and the pairing states. From analyzing the effect of the 
glide reflection, we are also able to show that treating the local 
spin-orbit coupling would require working with the 2-Fe FBZ. 
We should stress that we will consider the effect of the glide 
reflection in the context of the entire space group symmetry. 

The remainder of the paper is organized as follows. In Sec¬ 
tion II A, we aim to give a rigorous analysis of the mapping 
from the physical 2-Fe BZ to an effective 1-Fe BZ for a class 
of materials with P4/nmm space group symmetry. Our con¬ 
siderations apply to the ideally 2D case for which the con¬ 
duction electrons do not disperse along the c-axis. We show 
that in this case, the existence of a glide reflection symmetry 
for all momenta is guaranteed by this particular space group. 
Without a spin-orbit coupling, the classification of all Bloch 
states under the glide reflection and the particular degenera¬ 
cies along the BZ ensure that the folding procedure does not 
violate any space group symmetries of the system. In Sec¬ 
tion IIB we discuss some consequences of the space group 
symmetry on the superconducting pairing in the multi-orbital 
t, — Ji — J 2 model. In Section III, we discuss the effects of 
an atomic spin-orbit coupling term on the electronic proper¬ 
ties of both normal and superconducting states. We first show 
in Section III A that despite the glide reflection still being 
a valid symmetry for all momenta, the lack of degeneracies 
along part of the BZ edge nullifies the usual unfolding proce¬ 
dure. In Sec. IIIB we present and compare numerical results 
for the normal state bandstructure with and without a spin- 
orbit coupling to confirm the preceding symmetry-based ar¬ 
guments. The direct consequences of the spin-orbit coupling 
on the pairing are then discussed in Sec. IIIC. In Sec. IV we 
show how the symmetry arguments on the BZ folding survive 
for a 2D dispersion even when the interlayer couplings are 
taken into account. We also discuss the constraint imposed by 


the space group symmetry on the glide-reflection parity of the 
pairing channels, and subsequently examine the compatibil¬ 
ity of several proposed pairings with the previously described 
space group symmetry. Concluding remarks are given in Sec. 
V. Appendix A contains a derivation of the important ansatz 
used in the main sections while related details are given in 
Appendix B. 


II. THE P4/nmm SPACE GROUP SYMMETRY AND THE 
BRILLOUIN ZONE FOLDING 

In a large group of Fe-based materials which have 
the PA/nmm space group symmetry, the identical Fe- 
pnictogen/chalcogen layers are stacked on top of each other 
along the c-axis. The nontrivial spatial symmetry properties 
can be traced back to the structure of a single layer which is 
composed of a square Fe lattice in between two square As lat¬ 
tices shifted horizontally w.r.t. the Fe lattice and each other. 
The projection of the layer onto the Fe plane is shown in Fig. 1 
(a). Two adjacent Fe sites have different nearest-neighbor As 
configurations defining the A and B sublattices. 

The space group of the 1111 Fe-based superconductors is 
P4/nmm. The latter is non-symmorphic such that under any 
choice of unit cell one cannot decompose the set of symme¬ 
try operations into a point subgroup and its coset made up 
of proper lattice translations. In particular, for the conven¬ 
tional 2-Fe unit cell choice [14], the crystal is invariant under 
a glide reflection symmetry {u z \^0} composed of a frac¬ 
tional unit cell translation T t = {E | | ^0} in units of the sub¬ 
lattice translation, followed by a reflection about the Fe-plane 
P z = {a z \ 000}. The notation was chosen to be consistent 
with Ref. 1. 

The crucial point to consider is that the glide reflection 
T t P~ can be used to classify Bloch states of arbitrary momen¬ 
tum k. As observed in Ref. 2, for a general k = ( k x . k y . 0) 
in the Folded Brillouin Zone (FBZ) corresponding to the 2-Fe 
unit cell, the group of the wave-vector is isomorphic to the 
Cih point group. The latter has two irreducible representa¬ 
tions which have even/odd parity under the simple reflection 
P z . Consequently, states belonging to the irreducible repre¬ 
sentations of the space group P4/nmm for general k are also 
either Even (E) or Odd (O) under the glide reflection. In our 
view, this provides a connection between the particular glide 
reflection symmetry as part of the space group of the Hamilto¬ 
nian and it’s use in classifying states of arbitrary 2D momen¬ 
tum k , an argument we feel lacks an explicit exposition in the 
literature. By general group theoretical arguments [15], states 
belonging to different irreducible representations cannot mix 
and thus the even/odd Bloch states do not hybridize. As we 
show in the following, this, together with the particular degen¬ 
eracies along the entire folded BZ edge (see Fig. 1) allows the 
reduction to the 1-Fe unit cell. 
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FIG. 1. (a) Projection onto the Fe plane of the crystal structure of 
1111 systems. The thick lines define the crystallographic 2-Fe unit 
cell while the dashed lines denote the effective 1-Fe unit cell. As 
throughout the text, the unit distance is defined by the NNN Fe trans¬ 
lation (same sublattice), (b) The Folded Brillouin Zone (FBZ-thick 
line) corresponding to the true 2-Fe unit cell and the unfolded BZ 
(UBZ-dashed line) of the effective 1-Fe unit cell (for definiteness, in 
the k z = 0 plane). The FBZ is defined by the translation in recip¬ 
rocal space by Q = (±27t, 0,0) or (0, ±2-7 t, 0), which correspond 
to (±7r, ±7r, 0) or (±7r, ±7r, 0) in the notation for the effective 1-Fe 
unit cell. Here F, X and M label the points in the FBZ, while A, S 
and Y mark the corresponding segments. 


A. The glide reflection symmetry and the Brillouin Zone 
folding 

The advantage of using an unfolded 1-Fe BZ consists in ef¬ 
fectively reducing the number of bands in the calculation. In 
practice this means introducing a set of trial wavefunctions 
defined on an UBZ, carrying on the calculation, then folding 
back to the physical FBZ. In doing this, one must check that 
the Bloch functions defined in the UBZ do not introduce addi¬ 
tional hybridization terms in the Hamiltonian, and that they do 
not violate any of the original symmetries of the space group. 
Here, we analyze the same procedure in reverse. We start with 
a set of trial wave functions which are in accord with a mini¬ 
mum space group symmetry, namely the glide reflection, then 
examine the unfolding process. 

As mentioned in the previous section, the glide-reflection 
can be used to classify the irreducible representations of the 
space group for general k in the folded BZ when there is no 
dispersion along the z-axis i.e. k z = 0. This ensures that 
eigenstates of this operator will be eigenstates of the Hamilto¬ 
nian within a transformation on the orbital indices alone. As 
a consequence, the Hamiltonian can be written in block di¬ 
agonal form at these general k. The detailed representation 
theory of this space group shows that for a number of higher- 
symmetry momenta in the FBZ, the irreducible representa¬ 
tions are not required to be eigenstates of the glide [2], This 
does not invalidate our arguments since by continuity any off- 
diagonal terms must vanish here as well. We will show that 
this implies a one-to-one correspondence between the eigen¬ 
states of the glide-reflection and Bloch states defined on the 
UBZ, with no hybridization between different momenta. 

We proceed to build Bloch eigenstates of the glide opera¬ 


tion with k in the folded zone (See Appendix A). The elec¬ 
tron annihilation operator in the physical 2-Fe BZ on sublat¬ 
tice A(B) is defined as 

Ck,A/B,a = ~^^'Yl elk ' Ri ’ A/BG Ri,A/B,on ( 1 ) 

where a is an orbital index, N s = N/2 refers to the number 
of (2-Fe) unit cells and i is the index of the unit cell. The 
two sublattice vectors are related by Ri a + r = R,.n with 
t the nearest-neighbour distance used in the translation of the 
glide operation. The anihilation operators do not have definite 
parity under the glide reflection T t P z : 

C T T P z )C k , A/B , a = (-1 (2) 

where (—1)“ = ±1 depending on the parity of the local or¬ 
bital under a pure reflection. Without loss of generality we 
can define operators with definite parity : 

Cfe.B.a = [Ck,A,a + ( —1 )“<?*,B,a] (3) 

Cfc.O.a [^ fc A,a — ( —1 ) a Ck,B,a] ■ (4) 

Note that the even ( E ) and odd ( O ) parity states refer to parity 
under the glide and so we can build both types of operators for 
arbitrary orbital parity a. We also refer the reader to Appendix 
A for more details on the states defined above. 

Since Ck,E,a and (-'k.o.<> have different parity under the 
glide for arbitrary ”2-Fe” crystal momentum , they do not mix 
in a one-particle Hamiltonian. An arbitrary 2D tight-binding 
Hamiltonian consistent with the space group symmetries can 
then be expressed in the terms of the E/O states as 

Htb= Y, ('eW C 1 E, a Ck,E,p + e a o P (k)Clo ia C k ,o,p 

kGFBZ 

(5) 

where FBZ stands for the 2-Fe Folded Brillouin Zone and 
e if (k), tQ (k) are matrices in orbital space. We omitted the 
spin index for simplicity. 

We now turn to the unfolding procedure and start by noting 
that the eigenstates of the glide reflection have (see Eqs. A10 
and All) 

Ck+Q,E/0, a = Ck,0/E,a, ( 6 ) 

where Q = (±27t, 0) or Q = (0, ±27t) in units of the 2-Fe 
unit cell, and correspond to (±7r, ±7 r, 0) or (±7r, ±7r, 0) in 
units of the effective 1-Fe unit cell. That is, our states ap¬ 
parently violate the symmetry under the pure translation by a 
2-Fe unit cell sublattice. In the simplest Bravais lattice case a 
state of arbitrary momentum is only labeled by the irreducible 
representation of the pure translation subgroup i.e by k itself. 
Here, in addition to k one can also label all the states accord¬ 
ing to the irreducible representations of the glide reflection as 
well, that is, by both k and A. The resolution to our apparent 
conundrum lies in the fact that states displaced by a reciprocal 





lattice vector Q can switch their representation with respect to r = (±|,±|) in units of the 2-Fe unit cell. Explicitly, the 
the glide reflections: eigenstates of the glide become 


E^O,ef(k + Q)=eg(k). (7) 

This ensures that we are not violating the original pure transla¬ 
tion symmetry but it also introduces a constraint: Eigenstates 
of the system must be at least doubly degenerate all along the 
edge of the first (folded) BZ. Only in such case can the states 
switch representations as one goes beyond the folded BZ. In¬ 
deed, this is precisely what happens in Si and Ge at the X point 
along the edge of the BZ [16], [17]. The lattice in this case has 
a diamond structure and the two sublattices can be connected 
by a glide reflection. Due to the non-symmorphic nature of 
the space group, the states are degenerate at X and must ex¬ 
change representations as one goes into the extended BZ. In 
the present case for PA/nmm, there is a minimal degeneracy 
throughout the edge of the folded BZ [2], [18] which allows 
for the switching of the representations as one crosses into the 
extended zone. 

We can connect with one of the initial folding proce¬ 
dures [ 1 ] which relied on a definition of a quasi-crystal- 
momentum dependent on the parity of the orbital state under 
the pure reflection. Using the equivalent form of the glide 
states in (A10), (All), we can absorb the negative sign from 
(_ 1 )(« =1 ) = e ±iQ-r _ for states with odd parity under 
a pure reflection. Here Q = (±27r, 0) or Q = (0, ±27r) and 


Ck,E,o 


Cfc.a — 


1 

VN 




( 8 ) 


Ck,E,a odd = Ck+Q,a = ^ 7 ^ Y. e »(fc+Q)-fli (j R 


Ck 


,0 ,ck even Ck-\-Q ,ot 


1 £[ e < c +«>«, Clii0 

Ri 


(9) 


GO) 


C k 


,0,at odd 


Ck.a = ~r= 


1 

Yn 


Y [e ik ' Ri C Rz ,c 


( 11 ) 


The expressions above can be formally subtituted into the 
Hamiltonian (5). Under this mapping, the initial summation 
over k £ FBZ for both E and O glide parity sectors can be 
re-written as a summation over k £ UBZ (unfolded BZ) for 
the E only sector (recall that e^ /3 (k + Q) = e°p(k)). By 
taking into account the fact that different orbital parity states 
must be displaced by Q w.r.t. to each other we can recast the 
tight-binding Hamiltonian as 


Htb 


= £ £ 

leUBZ L ee 


+ £ f 'T't ^ t lQ l / f Q .) 


= £ £ ■ 

keUBZ «/5 


aCf+Q,#) + Y ( e a/3 ^) C l+Q, a C f: 

eo oe 

( 12 ) 

(13) 


where e, o refer to the parity of the orbitals a and ,3. The de¬ 
tails of the above derivation are presented in Appendix B. In 
going to the second line we explicitly identified along with 
Ref. 1 the dependence of k on the orbital parity such that 
k — k f° r even orbital states and k = ft + Q for odd states. 
What allows us to equate the two forms is the fact that ft is 
effectively a dummy variable in ( 12 ) and the matrix structure 
is determined by the orbital parities alone. For the purpose of 
notational clarity, we stress the following: a) k is defined in 
the larger Brillouin zone (UBZ) for the 1-Fe unit cell, and is 
associated with an orbital-dependent procedure in unfolding 
from the smaller Brillouin zone (FBZ). It has the same mean¬ 
ing as that used in Ref. 1; b) A;, also defined in the UBZ for 
the 1-Fe unit cell, arises in an orbital-independent unfolding 
procedure. The distinction between these two wave vectors 
will be important for our later considerations. 

We can illustrate a different widely used unfolding pro¬ 
cedure if we define two types of local states. Specifically, 


we can set Cr a . c , = Cr a i T .a f° r even parity orbital a and 
Cr a ,c = (—1 )CR A + Tta for odd parity orbital a in Eqs. 8-11 
such that the local orbital state changes sign under the sim¬ 
ple 1-Fe unit cell translation T t . The additional minus sign 
for odd orbital states will generate and additional e l ® T factor 
for odd orbital states. This is equivalent to the orbital parity 
dependent gauge transformation encountered in the literature. 
Using the same notation as above the eigenstates become 

Ck,E,a =Ck,a (14) 

Ck,0,a =C'k+Q,a (15) 

for arbitrary parity a. If we identify + Q) 

and use the above definitions in the Hamiltonian (5) the latter 
can be re-written as 
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E (4 /, (*)c£ 1< A,/s 


kGFBZ ' 


+ C h ' '.k + Q)(’1 + q a Ck+Q : /j 

(16) 

E 

(17) 


keUBZ 


where k is in the Unfolded Brillouin Zone (UBZ). 

We stress that both unfolding procedures discussed above 
were guaranteed by the absence of mixing between the E, O 
sectors in the Hamiltonian for general k together with the 
switching of the glide parity for the ansatz states beyond the 
first folded BZ (Eq. 6 ). These two conditions are a direct con¬ 
sequence of the irreducible representations of the glide sym¬ 
metry for arbitrary k G FBZ. None of the folding procedures 
discussed above and used in the literature would work if the 
Hamiltonian contained E, O mixing terms since we could not 
define Q shifted states related to the 1-Fe unit cell Fourier 
transforms for either even and odd parity orbitals. 

Although the above arguments seem to be implicit in the 
literature [ 1 , 12 ], we have encountered few observations of 
the crucial connection to the irreducible representations of the 
PA jnrnm space group, with a notable exception of Ref. 13 in 
the context of the electronic structure. Here, by constructing 
the electronic states which are even or odd under glide reflec¬ 
tion, we are in position to discuss the the constraints imposed 
by the glide-reflection symmetry on the various types of pair¬ 
ing, which will be discussed in the next subsection, as well 
as the effects of a local spin-orbit coupling on the unfolding 
procedure, which we turn to in the following section. As we 
mentioned in Sec. II, we believe that it is insufficient to argue 
that the presence of a glide reflection operation in the space 
group of the Hamiltonian guarantees the success of a folding 
procedure. Along with the glide reflection, the space group 
contains fifteen other symmetry elements [14] different from 
pure sublattice translations. Of these, only the glide reflection 
can be used to classify states according to the irreducible rep¬ 
resentations of the space group for general two-dimensional 
crystal momenta, since only this operation leaves an arbitrary 
k invariant. But as we will show in Sec. Ill, the same fold¬ 
ing procedure fails when a spin-orbit coupling is turned on, 
even though the additional interactions do not violate the glide 
reflection symmetry or any space symmetry of the original 
Hamiltonian. 


B. The effects of the glide symmetry on the superconducting 
pairing interactions 

We now discuss some of the effects of the glide symme¬ 
try for the superconducting pairing interaction. We assume 
that the symmetry outlined in the previous sections is not vi¬ 
olated such that without the pairing interactions one can un¬ 
fold the BZ of the 2-Fe unit cell to the effective 1-Fe UBZ. 
We also assume that the original symmetry of the 2-Fe unit 
cell is not broken by the appearance of magnetic or “ne¬ 
matic” order. The analysis presented here also assumes a triv¬ 


ial spatial dependence of the bare pairing interactions such as 
would arise in a nearest neighbor (NN), next-nearest neighbor 
(NNN) t — J\ — J 2 Hamiltonian. Lastly, we consider singlet 
pairing, as evidenced by experiments for various iron-based 
superconductors. In this case, the antisymmetric nature of the 
pairing wavefunction requires the pairing to be even parity un¬ 
der inversion [23], 

Before proceeding in a manner analogous to that of Sec. 
II A, some additional remarks are in order. The tight-binding 
part of a Hamiltonian with pairing interactions at mean-field 
level is typically chosen as the identity representation of the 
space group. In general, the pairing functions are determined 
self-consistently and can lower the symmetry of the Hamil¬ 
tonian to a subgroup of the full space group (as e.g., what 
happens to the rotational invariance of the Hamiltonian under 
the C 4 operation of the £> 4 /, point group of the full P4 /nmm 
space group in the case of a d-wave pairing). The problem 
simplifies for the spin singlet pairing considered in this sec¬ 
tion. For Cooper pairs of equal and opposite momenta in the 
FBZ ( k, —k ) associated with the 2-Fe unit cell, the symme¬ 
try properties of the pairing functions are completely deter¬ 
mined by the D±h point group associated with the PA/nmm 
space group. More precisely, the tensor irreducible represen¬ 
tations of the space group used to classify the pairing are the 
irreducible representations of the group of the wave-vector la¬ 
beled by total momentum k + (— k) = (0,0) (For a general ar¬ 
gument, see Ref. 22). Given this, we know that the irreducible 
representations of the D 41 , point group are either even or odd 
under inversion. In addition, for zero total momentum repre¬ 
sentations, the glide reflection is equivalent to a simple reflec¬ 
tion about the z-plane. This means that the parity under in¬ 
version coincides with the parity under the glide-reflection for 
one dimensional irreducible tensor representations and is op¬ 
posite for the two-dimensional ones [2], Restricting ourselves 
to the former ID representations, a Hamiltonian containing 
any linear combination of inversion-even pairing terms cannot 
break the glide reflection symmetry, regardless of whether the 
pairing preserves the rotational invariance or not. The glide 
symmetry can be broken if and only if the system sponta¬ 
neously breaks the inversion symmetry, which consequently 
results in triplet pairing. The remaining part of the section 
seeks to illustrate that no finite momentum pairing terms can 
be present in either the folded or unfolded BZ provided that 
the Hamiltonian is invariant under the glide-reflection. This 
also guarantees that the common folding procedure is still 
valid in this case. 

We consider for illustration purposes the following NN 
pairing interaction: 

Hint = £ £ J ap [cbt^m - c L c m x 

a/3 ij,NN 

~ Cia^Cj^ + H.C. 

(18) 

Based on Eqs. 3, 4 we can define Fourier transforms on 
each sublattice as 
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C RA , a = ‘ £ e -' fc ' Ra [C k , E , a + C k , 0 , a ], (19) 

V s keFBZ 

C RB , a = (-l) a -t= J2 e~ tkRB [Ck,E,a-C k ,o, a }. 

V s kGFBZ 

( 20 ) 


It is straightforward to see that the pairing interaction in Eq.18 
will contain equal numbers of E and O states. Indeed, this 
can be understood from the general discussion in the previ¬ 
ous sections. Although formally at some special points of the 
FBZ the eigenstates are not required to also be eigenstates of 
the glide reflection, in practice, continuity must enforce this 
as was mentioned in the preceding section. All terms in the 
Hamiltonian must effectively be invariant under the glide re¬ 
flection symmetry for arbitrary k , ensuring that all interactions 
must be even under the glide reflection symmetry. This con¬ 
strains all pairing terms to have equal numbers of even and 
odd states by the space group symmetry arguments. The gen¬ 


eral form of the interactions can in principle generate pairing 
terms mixing E and O states which might result in a finite Q 
momentum Cooper pairs. 

To make progress, we can consider (18) in a mean-field 
(MF) approach. In view of the connection between the MF 
Hamiltonian and the equation of motion approach we antici¬ 
pate the same results beyond the simplest level. As an illus¬ 
tration we analyze the term 


H Lt,MF - y. 

a/3 <ij> e 


(Jap [cU& m - Cl^ m ] ) X C^C m . (21) 


The remaining terms in Eq.18 can be obtained by flipping the 
spin indices of the last pair and adding the Hermitian conju¬ 
gate terms. Decomposing the NN summation as '}2 <i j > = 

EiU E e=±x,y + E Rb E e=±x,y and takin g the sublattice 
specific Fourier transformation (F.T.) as in (19),(20) we ob¬ 
tain 


H'int,MF — e ke x 

a/3 e keFBZ 

(—l) /j Ag 4 a/3 (Cfc i £) )a |C , _fc i £;,/34. ~ Cfc,-E,atC-fc,0,/34. + Cfc.O.a fC-k,E,0l ~ Ck,0,atC-k,0,/3l) + 
+ (—l) a Ae B , a p(Ck,E,at < ~'-k,E,l3l + Ck,E,oi^C~hp\.0 — Cfc,0,atC-k,.E,/3.t ~ Ck,0,atC-k.O,Pl) 


( 22 ) 


where 


A 


A 

e,a/3 


A 


B 

e,a/3 



° Ra at Ra +e/3i ° R a a±^R A +e0t 

(23) 



Rb& t -RB+e/3-t Rj3a\. i?s+e^t 


(24) 


The crucial assumption made in the above equations is that 
both A f, Af are independent of their sublattice space indices 
Ra and Rr , but A) 1 / Af in general. This amounts to 
having the pairing order parameters which are constant on the 
respective sublattices. 

We distinguish two cases: i) a and /3 corresponding to or¬ 
bitals of the same parity under the reflection, ii) a and /3 have 
different parities. 

For case i) of same parity we can look at the real space ex¬ 
pression (21) and demand invariance under a glide reflection. 
Since both a, /? terms acquire the same orbital parity factor 
and each state gets shifted by one unit of the 1-Fe unit cell 
we have A^ a p = A^ Q/3 . Plugging this into (22) and setting 
(—1)“ = (—1)^, we see that all EO and OE mixed terms 
cancel and we get only same glide parity terms. 


For case ii) where a and /3 have different orbital parities 
we can do the same as above and impose the glide reflection 
symmetry in real space. We get A^ a p = —A// a p. With 
(—1)“ = —(—l)' 3 we get the same cancellation as for case i). 


Identical results are obtained for the next-nearest coupling 
(NNN, J 2 ) case. 


In both cases of same and different orbital parity pairing, 
the results are completely analogous to those considered in 
Section IIA for a tight-binding model without spin-orbit cou¬ 
pling. Namely, the pairing part of the Hamiltonian splits into 
EE and OO glide parity sectors. The exact same unfolding 
arguments can be trivially extended to the pairing part. Calcu¬ 
lations can be done in an UBZ using the k F.T. and then fold 
the results to the 2-Fe BZ. In this scheme there is no finite 
momentum pairing. 
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III. THE EFFECTS OF SPIN-ORBIT COUPLING 
A. Folding in the presence of spin-orbit coupling 

A local (atomic) spin-orbit coupling (SOC) term 

IL- S = I -(J 2 ~ L 2 ~ S 2 ) (25) 

preserves the PA/nmm space group symmetry since the lat¬ 
ter is a scalar under all point group operations and is local 
in space (I is a constant). It does however lock the orbital 
and spin parts of the conduction electrons together such that 
the two do not transform independently of each other under 
space group operations. This forces us to change the irre¬ 
ducible representations of the space group to double valued 
representations [18]. In general, the reflection P z = {<r z |000} 
is equivalent to 

a z = C 2z ® I (26) 

where C-> z corresponds to a rotation by n along a chosen z- 
axis and I is the inversion. Since spinor states are invariant 
under the inversion [18], the effect of P, is given by the Pauli 
matrix term 



(27) 

For arbitrary momentum k , one constructs the irreducible 
representations of the space group by determining the so- 
called group of the wave vector. That is the subgroup of the 
space group, with elements which either leave k invariant or 
translate it by a reciprocal vector. By examining the effect 
of these elements on a set of properly chosen states, one can 
determine the group isomorphic to the group of the wave vec¬ 
tor. In the case of double valued representations one must 
also consider, in addition to the operations in the single-valued 
case, those obtained by changing in sign of the states [18]. 

Even with SOC terms, the glide reflection is still part of the 
group of the wave-vector. Because of the locking of the spin 
and orbital states, the glide reflection generates factors of i, 
consistent with Eq. 27. By considering local states with and 
without spin and applying the operations that keep an arbi¬ 
trary k invariant one obtains a group isomorphic to the double 
valued group of Cih- The irreducible representations of the 
latter are illustrated in Table III A. 

The inclusion of the spin in the glide reflection will always 
generate the pure phases ±i. We then conclude that the phys¬ 
ical irreducible representations correspond to T 3 and T^ This 
indicates that, as for the PA/nmm case without SOC, for ar¬ 
bitrary momentum in an unfolded BZ, eigenstates of a general 
tight-binding Hamiltonian will also be eigenstates of the glide 
reflection symmetry. We can connect with the ansatz states in 
Eqs. 3, 4 which coincide with Bloch states in a 1-Fe BZ. We 
remark that these are still eigenstates of the glide reflection 
operator but in the presence of SOC, we need to account for 
the transformation of the spins as well. We thus relabel 


TABLE I. The double-valued irreducible representations of the Cih 
point group [ C.J. Bradley and A.P. Cracknell, Mathematical Theory 
of Symmetry in Solids (Clarendon Press, Oxford 1972) ]. 

Cih E o z Ed z 

IT 1 1 1 X 

r 2 l -l i -l 

IT 1 * —1 — i 

IT 1 —i —1 i 


Ck,E,ct ,t Cfc.O.a,t (28) 

Ck,E,a,l ~> ^k,E,a,i (29) 

Ck,Q,a ,t C k E ^ t (30) 

Ck,0,a,i -A Ck,6,a,i ( 31 ) 


where the E, O refer to the sign in front of the ie lk ' T term 
under the glide transformation. Since the irreducble repre¬ 
sentations are also eigenstates of the glide, the Hamiltonian 
excluding the pairing terms can only connect E, E or 0,0 
states. Note however that states E , O states in the original (no 
SOC) labeling such as C'k.E,a ,fj Cfc.o,a,t b°th belong to the 
E irreducible representation of the space group with SOC. So 
in general, the space group symmetry allows the mixture of 
the E, O states invalidating the unfolding procedure since the 
two will always correspond to states with shifted momenta k 
and k + Q. 

Also note that the eigenstates of the Hamiltonian in this 
case are not degenerate along the Y line of the unfolded BZ 
[2], [18] (not counting Kramers degeneracy which is irrele¬ 
vant here). In general, there cannot be an analogous switching 
between the T 3 and T 4 representations at the edge of the BZ 
and one must conserve the “parity” under the glide reflection 
as one crosses into adjacent zones. Therefore for finite SOC 
and for a general choice of tight-binding parameters, 

C'k,E/0,a — Cfc+Q,0/£,a , (32) 

with ±z eigenstates under the glide-reflection for T 3 and T 4 
respectively. This forces us to accept other eigenstates of the 
glide-reflection such as those in Eqs. A14, A15. 

The relation above also signals that we cannot choose the 
simple sublattice superposition states in (3),(4) and so we can¬ 
not connect with the effective extended momentum Fourier 
transforms in Eqs. 8-11 or the gauge transformed states in 
(14),(15). In a practical calculation where the 1-Fe unit cell 
assumption is of any use we implicitly carry on calculations 
using the familiar 1-Fe unit Fourier transform defined on 
k £ UBZ and then fold to the physical 2-Fe unit BZ. But 
our arguments show that in doing so we are violating space 
group symmetry and as a consequence we have no guarantee 
that the results thus obtained coincide with those done directly 
in the unfolded zone. 

As in the zero SOC case, we argue that it is insufficient to 
justify the validity or invalidity of the 1-Fe unit cell by invok¬ 
ing a general symmetry of the Hamiltonian, namely, the glide 
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reflection. In the SOC case, the space group of the Hamilto¬ 
nian does not change since the extra coupling does not violate 
any symmetry of the latter. Our arguments, through which 
we attempt to treat the glide reflection within its natural space 
group symmetry perspective, can be used to give a more pre¬ 
cise prescription to its use and to illustrate how the formula¬ 
tion of the problem in an 1-Fe unit cell can fail in spite of the 
validity of the glide reflection as a symmetry of the Hamilto¬ 
nian. 


B. The effects of the spin-orbit coupling on the normal-state 
bandstructure 



Uk_!_L_I_i_!_I_ 

r k x /jt M 


FIG. 2. (a) Bandstructure of the five-orbital tight-binding model in 
the folded Brillouin zone (FBZ) without the spin-orbit coupling, (b) 
The corresponding Fermi surface in two quadrants of the FBZ. 

In light of the discussion on the folding procedure, in this 
section we study the effects of the spin-orbit coupling on the 
bandstructure of the normal state. We consider the following 
Hamiltonian: H = Htb + Ffscn where 

Ht b = E ^(k )Ci aa C^ ai (33) 

ka/3,£7 

is a five-orbital tight-binding model for the parent compound 
of iron pnictides. Here e ct) g(k) are tight-binding parameters, 
which we adopted fromRef. 19. //so = Aso Xo k,-S,, refers 




FIG. 3. (a) Bandstructure of the five-orbital tight-binding model 
in the folded Brillouin zone (FBZ) with a local spin-orbit coupling 
Aso = 0.05 eV. (b) The corresponding Fermi surface in two quad¬ 
rants of the FBZ. 


to a local spin-orbit coupling term. As we discussed, this term 
does not break any symmetry of the lattice, but couples the 
spatial and spin part of the single-particle wave function. 

According to Sec. II A, the tight-binding Hamiltonian de¬ 
fined by the pseudo-crystal momentum k in the UBZ can be 
transformed to the physical momentum in FBZ via the folding 
procedure. This allows us to study the effects of the spin-orbit 
coupling by comparing the bandstructures without and with 
a spin-orbit coupling in the FBZ. In the absence of the spin- 
orbit coupling, the bandstructure of the tight-binding model 
is shown in panels (a) and (b) of Fig. 2; it is in agreement 
with that from ab initio calculations using the density func¬ 
tional theory [19]. It is clearly seen that the bands are doubly 
degenerate along the boundary of the FBZ (Y line from the 
M to the X point), and as a consequence, the two elliptical 
electron pockets cross at a point along this direction. As we 
have emphasized in Sec. II A, this double degeneracy guar¬ 
antees the successful folding procedure: the wave function 
can switch representations under the glide when crossing the 
FBZ boundary, and hence one can define a pseudo-crystal mo¬ 
mentum according to the parity of the wave function under 
the glide. This ensures the equivalence between the models 
defined in the FBZ and UBZ. When a spin-orbit coupling is 
turned on, as shown in Fig. 3 (a) and (b), the double degener- 
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acy along the FBZ boundary (M-X direction) is lifted. As we 
discussed, this invalidates the unfolding procedure since the 
wave functions have to conserve the parity under glide across 
the FBZ boundary. As a result of the lifted degeneracy, the 
two electron pockets no longer cross, but a gap opens along 
the X-M direction. 

Another feature of the bandstructure in the presence of spin- 
orbit coupling is the opening of hybridization gaps between 
the d xy and d xz / yz bands along the 1-M direction as shown 
within the dashed ellipses in Fig. 3. In absence of spin-orbit 
coupling, these bands simply cross, without opening a gap 
because they have different parity (pseudo-crystal momenta) 
and can not mix. But with a finite spin-orbit coupling, the 
C^ xz/yz,^ an d fk+Q,£cy,| have the same parity and they can 
hybridize via opening a gap. Recently, a hybridization gap 
between the d xz / yz and d xy bands along the T-M direction has 
been observed in ARPES measurements on several iron-based 
compounds [20, 21], It would be interesting to compare the 
experimental results with theoretical ones as this may provide 
valuable information about the strength of spin-orbit coupling 
in these systems. The spin-orbit couplin g also mixes the d xz 
and d yz orbitals, and lifts the double degeneracy between the 
d xz and d yz orbitals at V point. But a local spin-orbit coupling 
does not lift the degeneracy between the d xz and d xy orbitals 
at the bottom of the electron bands at M point since it does not 
break the four-fold rotational symmetry of the P4/nmm group. 


C. The pairing interactions in the presence of spin-orbit 
coupling 

In Sec. Ill A we showed that the correspondence between 
states defined on a 1-Fe BZ and those belonging to the irre¬ 
ducible representations of the space group in the presence of a 
SOC breaks down even when there are no pairing interactions 
present. For consistency, here we analyze the direct effect of 
a local spin-orbit coupling term (Eq. 25) on the pairing inter¬ 
actions of a t — J model. 

The projective nature of this model [3] excludes non-singlet 
pairing terms. However, the SOC does not conserve the elec¬ 
tronic spin quantum number and in this case a t — J model 
must include ad-hoc triplet pairing terms. In general, the re¬ 
sulting pairing must include terms which are odd under inver¬ 
sion. Based on the discussion at the beginning of Sec. IIB, 
we expect that the inversion-odd pairing, in the spin triplet 
case, will also be odd under the glide-reflection and thus cor¬ 
respond to finite-momentum pairs in the UBZ. We note that, 
as in the case for pairing without SOC, we restrict ourselves to 
one-dimensional representations of the point group, for which 
the above is correct. In the case of two-dimensional repre¬ 
sentations, the opposite holds, with inversion-even functions 
actually corresponding to glide-reflection odd states. 

At mean-field level, the pairing can still be written in the 
form of Eq. 21 except for the spin structure allowing both 
same-spin and opposite spin pairing. More specifically, we 
use the correspondence in Eqs. 28-31 between the E, O la¬ 
bels without SOC and the E, O ones with SOC turned on to 
re-write the pairing in Eq. 22 in terms of the latter. As a 


consequence, the real-space pairing functions A^ 1 ^, A R a/3 
(Eqs. 23, 24) acquire a spin-index dependence. The previous 
argument (without SOC) relied on the transformation prop¬ 
erties of the latter term under the glide reflection. We can 
apply the same procedure in real-space allowing the transfor¬ 
mation of the spins. For singlet pairing the f, I combination 
always generate products of i, —i, ensuring that A a ’ (3 have 
exactly the same properties as before. Applying the transfor¬ 
mation (28)-(31) in reverse, we note that singlet pairing can 
only have diagonal E, E and (). O pairing (original labeling) 
and so does not introduce any finite momentum Cooper pairs. 

The situation is different when we allow triplet pairing. 
The same-spin terms always generate a minus sign under the 
glide reflection. When a, (3 have the same parity this means 
A = — A^ q/3 . In the original labeling, the only terms that 
survive are the off-diagonal E, O. Similarly, when a, (3 have 
different parity A^ ^ = A R a/3 but the (—1) terms coming 
from the orbital parity guarantee that we again obtain only off 
diagonal E, O pairs. Now however, the E, O terms will map 
to k,k + Q in the UBZ producing finite-momentum Cooper 
pairs and invalidating the unfolding. The spin-symmetric 
combinations allow forms like those for singlet pairing and 
thus do not introduce any new terms. 

IV. DISCUSSIONS 

In this section we elaborate on some further issues related to 
the space group symmetry and BZ folding we detailed above. 

A. Effects of three-dimensionality vs. spin-orbit coupling 

Real materials with a PA/nmm space group have a 3D 
structure ensuring the bands are always dispersive, albeit 
weakly, along the k z direction. The group of the wave- 
vector for arbitrary k z cannot contain the glide reflection since 
the latter connects the generally inequivalent k z , — k z compo¬ 
nents. For k z = 0, the folding can still work as we illustrate 
below. The dispersion along the ^-direction can be accounted 
for by generalizing the 2D ansatze (3), (4) to 

Ck,k z ,E,a =— / ^'^2 elkzRz \Pk,A,ot,R z + ( —1 ) a Ck t B,a,R z ] i 

(34) 

Ck,k z ,0,a = “/= e lkzRz [Ck,A,oc,R z — ( —1 ) a Ck,B,a,R z ] i 

v 2 R z 

(35) 

which for general k z are not eigenstates of T t P z and we 
keep the E, O indices for labeling purposes. Here, k still 
refers to a purely 2D wave vector. Since for finite k z , these 
states do not necessarily coincide with the irreducible rep¬ 
resentations of the space group, the Hamiltonian will con¬ 
tain terms mixing E, O indices. For k z = 0 however, the 
glide does not affect the Bloch momentum label but it maps 
R z —> —R z . Due to the equivalence of the different planes. 
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Ck,A/B,a,Rz = Ck,A/B,a.-R z since the two Wannier states 
are related by a proper translation. This means in the k z = 0 
plane the states (34), (35) still coincide with the irreducible 
representations of the space group and the previous arguments 
still apply. Therefore, states with two-dimensional crystal mo¬ 
menta can still be mapped onto an effective 1-Fe BZ even 
when the dispersion along the z-axis is turned on. 

B. Space group symmetry and parity of pairing 

Recent studies have paid particular attention to the conse¬ 
quences specific to the 2-Fe unit cell such as the use of an ad¬ 
ditional isospin quantum number [6] on the superconducting 
pairing. But some important issues are still unclear. For ex¬ 
ample, pairing channels with very different symmetries have 
been proposed theoretically, [6, 8] although the target com¬ 
pounds have the same spatial symmetry. As we discussed ear¬ 
lier in our paper, the space group already imposes strong con¬ 
straints on the parity of symmetry-compatible pairing. This 
provides a means to check whether a pairing channel is al¬ 
lowed by the space group symmetry of the system. Here we 
examine this issue for two recently proposed pairing chan¬ 
nels, the TAO pairing [6] and the //-pairing [8]. As before, we 
focus on spin-singlet pairings in one-dimensional representa¬ 
tions henceforth. 


1. TAO pairing 

Discussed in Ref. 6, the TAO pairing refers to a spin sin¬ 
glet, orbital triplet A\ g d x 2_ y 2 + id xy pairing in the 2-Fe BZ. 
It is easy to check that this pairing channel has an even parity 
and is compatible with the Pd/nmm space symmetry. From 
Sec. IIB we know that for a 2D dispersion of the conduc¬ 
tion electrons and in the absence of spin-orbit coupling, any 
pairing defined on the 2-Fe unit cell FBZ which respects the 
space group symmetry must have an equivalent, albeit given 
by a different linear combination of channels, (k, k) pairing 
in the 1-Fe unit cell UBZ. The aforementioned TAO pairing 
defined in the 2-Fe BZ (k, —k) is equivalent to one defined in 
the 1-Fe UBZ and does not incorporate any particular prop¬ 
erties of the 2-Fe BZ not captured by the former. A possible 
advantage of the direct 2-Fe unit cell formulation might con¬ 
sist in expounding physical features which might be harder to 
illustrate in an equivalent 1-Fe unit cell picture. 

2 . rj-pairing 

Recent discussions have also considered the //-pairing in 
the iron-based superconductors [8-11]. In its original pro¬ 
posal, [8] the //-pairing refers to a singlet pairing of two elec¬ 
trons with pseudo-crystal momenta k and k + Q. respec¬ 
tively. This pairing has nonzero total pseudo-crystal momen¬ 
tum in the 1-Fe UBZ, and the momentum dependent part of 
the wavefunction has odd parity under inversion. This pro¬ 
posal was based on the observation that the inversion center in 


an Fe plane lies half-way in between two (inequivalent) sites. 
In a real space basis, the inversion operation interchanges the 
two positions on the different sublattices. In a simpler case 
where the two lattices are equivalent, the antisymmetry under 
exchange of the pairing wavefunction forces the spatial In¬ 
dependent) part to be even under inversion. The existence of 
the two inequivalent sites opens up the possibility of odd par¬ 
ity momentum dependence since the overall spatial part could 
be described as a direct product of a purely ^-dependent part 
(in the FBZ, 2-Fe unit cell description) and a pseudo-spin ma¬ 
trix which captures the effect of the different sublattices. The 
purely ^-dependent part can have odd parity under inversion 
as long as the remaining degrees of freedom compensate with 
a minus sign such that the total inversion parity is still even. 
Although not explicitly stated in Ref. 8, such a state could cor¬ 
respond to a E,0 pseudo-spin singlet. Upon unfolding, the 
combination will generate the finite momentum (fc, k + Q) 
pairing. As already shown in Sec. IIB, for states of definite 
parity under inversion, or equivalently, of definite parity un¬ 
der the glide-reflection, the pairing can be written in general 
as a linear combination of terms with and without E, O mix¬ 
ing together with their respective momentum-dependent parts. 
In addition, if we do not break the glide-reflection symmetry, 
the total pairing must be even. In a purely two dimensional 
BZ, the glide-reflection cannot change any of the momenta of 
the pair, such that terms with E, O mixing change sign under 
the operation while those without do not. Indeed, as detailed 
in our Sec. IIB , all mixing terms must vanish, ensuring that 
upon unfolding no finite momentum pairs can be generated. 
We stress that this conclusion holds whenever inversion is not 
spontaneously broken for a two dimensional BZ. 

In recent studies on the //-pairing, there also seems to be 
confusion in the definition for this type of pairing in the lit¬ 
erature which, we believe, is associated with different subse¬ 
quent definitions of the 1-Fe BZ. While in the original work 
the //-pairing referred to a (k, —k + Q) state, in more recent 
studies, it is associated with a (ft, f + Q) pairing in the so- 
called “physical extended BZ” [10, 11] given by a mapping 
similar to the one in Eqs. 8-12. As in the original proposal 
[8], Ref. 10 alludes to the inequivalence of the two sublattice 
sites but does not attempt to explicitly consider eigenstates of 
the glide reflection operation. Rather, the authors classify 2D 
states both in the "physical” (f in our convention) and un¬ 
folded (k) representations according to the parity under the 
reflection about the z-plane. In the physical representation, 
an alternating minus sign in the hopping between different or¬ 
bital parity Wannier states is ’’absorbed” into the definition of 
the quasi-crystal momentum, resulting in a shift by Q between 
the even and odd reflection parity states, as detailed in our Sec. 
II A. The authors argue that a zero mometum pair in the UBZ 
(k. ■ k) or arbitrary parity under the z-reflection alone must 
decompose into a linear combination of even/even, odd/odd 
and even/odd terms which correspond to (/ + Q. 1ft — Q), 
{If, —If) and (ft + Q- —f) respectively in the physical BZ. 
From this, they seem to explain inconsistencies in the spec¬ 
tral weights in the superconducting state between calculations 
done on a 1-Fe/UBZ (k) and subsequently folded down and 
ARPES experiments among others whose results are naturally 
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obtained in a 2-Fe/FBZ. 

In the context of a glide reflection symmetry, the two con¬ 
structions referring to a finite momentum (k, —k + Q ) in the 
UBZ [8] and in the ’’physical” (ft, ft + Q) [10, 11] are dif¬ 
ferent. The (k. k + Q) pairing has an odd parity under glide 
reflection, while the (ft, —ft 4- Q) -pairing is parity even un¬ 
der the glide. To see this, consider the latter (’’physical” 
representation) case which corresponds to an odd/even par¬ 
ity combinations under the pure reflection. In our language, 
this corresponds to states given in Eqs. 8-11 with a odd and 
even respectively. It is not difficult to see that upon con¬ 
verting to the ’’true” 2-Fe/FBZ (k) this term corresponds to 
(k. —k ), O/O under the glide-reflection pairs. To the best of 
our knowledge, this point has not been clarified before. Ac¬ 
cording to our analysis, the only symmetry-allowed pairing is 
the even glide parity one, i.e., the (k. k) pairing (or equiv¬ 
alently, the (ft,— ft + Q) pairing). It corresponds to normal 
zero-momentum pairing in both the 2-Fe FBZ (k) or the 1- 
Fe UBZ in pseudo-crystal momentum space (k). More pre¬ 
cisely, the three formulations alluded to above are equivalent 
and we see no reason why 1-Fe unit cell/UBZ calculations 
which are folded down could not a priori capture experimen¬ 
tal results. This correspondence naturally explains why the 
superconducting gap functions obtained from calculations in 
the 2-Fe FBZ are identical to the previous results in the 1-Fe 
UBZ [11], Within our approach, we have also shown that the 
odd glide parity ( k, - k + Q) pairing is not allowed by sym¬ 
metry. 


singlet pairings in one-dimensional representations). With 
this criterion, the //-pairing with (k. —k + Q), originally pro¬ 
posed in Ref. [8] is not symmetry-allowed since it is parity 
odd. The //-pairing discussed in most recent works [10, 11], 
on the other hand, refers to a (ft. —ft + Q) pairing, which cor¬ 
responds to a total momentum zero ( k , - k ) pairing with an 
even parity, and is thus compatible with the PA/nmm space 
group symmetry. 

Our analysis for the folding within a PA/nmm space group 
symmetry can serve as a comparison point for a similar dis¬ 
cussion in the more involved /4 /rnrnrn case. As we stressed 
throughout the text, the validity of the glide-symmetry, to¬ 
gether with the fortuitous double degeneracy along the edge 
of the FBZ guarantee the success of the folding for k z = 0 and 
no SOC. Since none of these appear to be valid for lA/mmm, 
the folding will probably not work. A rigorous analysis of this 
latter case is reserved for a future publication. 
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Y. CONCLUSIONS 


The glide reflection symmetry is valid for states of arbitrary 
momentum and without spin-orbit coupling there is a minimal 
double degeneracy all along the 2-Fe unit cell BZ edge. This 
ensures that a tight-binding Hamiltonian can be determined 
using an unfolded BZ corresponding to a 1-Fe unit cell. By 
contrast, although the glide symmetry still holds for arbitrary 
momentum when spin-orbit coupling is turned on, the latter 
mixes states corresponding to different pseudo-momenta in 
the unfolded BZ and lifts the degeneracy along the Y line, 
forbidding the general use of the same unfolding procedure. 
These conclusions are consistent with bandstructure calcula¬ 
tions with and without spin-orbit coupling which show the 
lack of and the presence of this hybridization. 

We also conclude that for at — J type Hamiltonian with¬ 
out spin-orbit coupling the results obtained directly from a 1 - 
Fe unit cell should coincide with those from 2-Fe unit cell, 
the two being related by the validity of the unfolding proce¬ 
dure. This applies to the TAO pairing [6]: Though proposed 
in the 2-Fe BZ, it is equivalent to a d x 2_ y 2 + id xy pairing with 
both intra- and inter-orbital contributions. One more remark is 
that this equivalence does not hold when a spin-orbit coupling 
term is included. 

Another conclusion from our symmetry analysis is that the 
pairing channel compatible to the PA/nmm space group sym¬ 
metry must have an even parity (once again.we focus on spin- 


Appendix A: Eigenstates of the glide operation 

The irreducible representations of the P^/ nmrn space group 
for general k £ FBZ and for the special loci F, A, E and M 
[2] transform as 


(T r P z )C ka = e lkT X C ka , (Al) 

where A = ±1. The above form is not the case for X and 
Y. However, since the Hamiltonian must evolve continously 
with k it is clear that we can form irreducible representations 
at the two above mentioned loci by taking linear combinations 
of even ( E ) only or odd only ( O ) such that our arguments are 
not affected. 

We can derive the general form of states which transform 
according to (Al) from a general superposition of operators 
defined on each sublattice 




e lk RA C ( r A J o + e i9k e ikRs C\ 

Ra Rb 


( B ) 

Rboc 




flk-RA 


Ra 


(A) , c i8 k ikTf~t(B) 


C RAa + c k (-C 


Ra +ra 


(A2) 


where Ra, R/; are summations over the position vectors of 
the respective sublattices, and a stands for the parity of the 
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local degrees of freedom under a pure reflection a z . We stress 
that, in the most general case, the completely local states can 
be different ( Cra^ 7 ^ C'ra b due to the physical inequivalence 
of the two sublattice sites. Although in principle the linear 


combinations of the two sublattice Bloch states can have arbi¬ 
trary complex coefficients, the form chosen above is sufficient 
for our purposes. 

Applying the glide reflection to the trial state (A2) and im¬ 
posing (A2) we get 


1 


{T r P z )C ka = Ae' —— > e 


E e ikR A 

. 


Ra 


(“I rCRA+r* 


+ e i 6 k e ikT (-l) a C ' RA+2Ta 


,(B) 


(A3) 


such that all position vectors get shifted by the fractional 
translation r, the local states generate the parity term a, and 
the phase 6 k —> 6 k in general. Comparing (A2) and (A3) we 
see there are a number of possibilities. In general, condition 

(A3) cannot determine all the unknowns i.e. 

the phase fac- 

tor and the relation between the the two displaced local states. 
We can connect with a folding procedure by making some as¬ 
sumptions regarding the phase and letting the above condition 

determine the local states. 


A first possibility corresponds to taking A = 

1 in (A3) 

/'■((A) _( -1 \a/^(A) 

L Raoc ~( L > L Ra+tcx 

(A4) 

_ j 1 \a/~<(B) 

L R A +Ta ~V L ) L R A +2ra 

(A5) 

> ^ 

II 

(A6) 

For A = — 1 we can have 


/^<(A) / 1 ■va/'AA) 

L Ra a ~ ( L > L Ra +ra 

(A7) 

/~<{B) _ / 1 \<xsa(B) 

L Ra+toc ~ \ L > L Ra+2tu 

(A8) 

> * 

II 

(A9) 


and we set 6 k = 0. The conditions above are simply the trans¬ 
formation properties of the local states under a simple 1-Fe 
unit cell translation. If we choose the states at corresponding 
to ('a and C/> to have the same functional form, we can con¬ 
struct two distinct eigenstates of the glide by for arbitrary a 


C k , E , a =—Y j e lkRA \c RAal + (-1 ) a e ikT C RA+Ta 


Ra 


(A10) 




D ik-R.A 


Ra 


Cr ao - (-l) Q e ifcT C' flA+Ta j 
(All) 


where we omitted the A,B superscripts which are now irrele¬ 
vant. Eqs. A10 and All are the operators in (3) and (4). 

Note that in addition to the above states which allowed the 
unfolding we can also choose eigenstates of the glide reflec¬ 
tion by imposing 


CrI = ± (-1 ) a e^e‘ kT C^ +2Ta (A12) 

=± (A13 > 
We can choose the phase factor e l ° k = e l0k = e~ tkT and 
disregard the A, B distinction as in the first choice. The two 
sublattices are now related by a ^-dependent phase rather than 
a simple sign. The initial ansatz (A2) becomes 


Ck,E,Oi 


+ (-1 ) a C RA+Ta ] (A 14) 


Ra 


Cfc,0,a 


-(-1 ) a C RA+Ta \ (A15) 


Ra 


which clearly do not correspond to 

a 1-Fe unit cell. 


Appendix B: The Hamiltonian in the physical extended 
momentum basis 


Direct substitution of the definitions (8)-(ll) into the 2-Fe 
BZ Hamiltonian (5) gives 
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htb — e 


keFBZ 


E k )°L C w + 0*) C jLq,« C *+<w) + E ( fc )^ + Q, a ^ +Ql/3 + e%(k)Cl a C k ,p 


+ E { e ap( k ) C l,a C k+Q,P + e ap( k ) C l+Q, a C k,pj + E ( e f/3( fc ) C Jfe+Q,a C '‘M + e ap( k ) C k, a Cfc+Q,/3 


(Bl) 


We can use the fact that efg 1 11 (k + Q) = e V> f (A) to rewrite We can also show that the third term can be expressed as a 

the first two terms as a total sum over Ijt over the unfolded BZ. sum over Specifically, 


E *Ut) C i aP C t + Q^= E ^ C l a ^ + Q^+ E tUt) C Lp C * + Q,<*0 (B2) 

IteUBZ IteFBZ I$fbz 

= E e a/3( fc )^la/3 C, fc+Q,a/3 + E 6 ap( k + Q) C l+Q } ocP C k+2Q, a p (B3) 


keFBZ 


keFBZ 


which after recognizing that C'k - iQ = Cfc gives the form in 
Eq. 12. 
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